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Abstract 

Contributions from high-order (non-perturbative triple and quadruple) cluster amplitudes to the 
dissociation energies, equilibrium distances, and vibrational constants for the ground states of van 
der Waals dimers Hg2 and Cn2 are evaluated. The incorporation of these contributions into the 
results of large-scale CCSD(T) calculations leads to non-negligible corrections of the computed 
molecular constants for Hg2 (6% for the dissociation energy), and enables one to attain perfect 
agreement with the experimental values. 
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I. INTRODUCTION 



The Hg dimer was the subject of a series of ab initio calculations (see Table [T]). It is 
well understood from these studies that only highly-correlated relativistic calculations can 
give reliable results for such systems as the Hg2 molecule in its ground state. The widely 
used DFT method with the popular B88P86 and PW91- exchange-correlation functionals 
cannot ensure an acceptable accuracy of the description of the Hg-Hg bonding^. In turn, 
within the wavefunction-based approaches, the basis set superposition error (ESSE) can 
seriously deteriorate the results for such weakly bound molecules even for rather extensive 
basis sets. To suppress the ESSE, large basis sets and counterpoise-like corrections should 
be employed. Fully relativistic calculations appear to be computationally too demanding 
for highly-correlated treatments and can be performed only with unacceptable restrictions 
for the basis sets^. The calculated spectroscopic properties of Hg2 reported in the papers^i^ 
are in a pretty good agreement with the experimental data. In Ref.-, the scalar relativis- 
tic coupled cluster method with single, double and non-iterative triple cluster amplitudes 
(CCSD(T)) and the (9s,8p,7d,3f)/[7s,6p,4d,3f] basis set were used to correlate the outermost 
24 electrons {5d and 6s shells). Then the effects of correlations with the 5s and 5p shells and 
the basis set extension to higher angular momentum functions were estimated at the MP2 
level of theory. Spin-orbit interactions were neglected at that stage. In Ref.-, 24-electron 
CCSD(T) calculations with an uncontracted (9s,9p,8d,4f,2g) basis set were performed. The 
correlation contributions from the 5s and 5p shells were estimated at the MP4 level of theory. 
Then the spin-orbit contribution was taken into account using MP2 calculations. CCSD(T) 
calculations in Ref.-- with the larger basis set (lls,10p,9d,4f,3g,2h) correlating 40 electrons 
(5s, 5p, 5d, and 6s shells) and with the spin-orbit correction taken from Ref.- are in a worse 
agreement with the experimental data. As has been noted in Ref.~, the good agreement of 
the calculated spectroscopic properties in papers^*^ with the experimental data is a result of 
the fortunate cancellation of errors. Our CCSD(T) calculations^ with the generalized corre- 
lation (14s,12p,9d,5f,3g,2h)/[10s,9p,7d,5f,3g,2h] basis set and spin-orbit contributions from 
the density-functional theory give results comparable with those in the work^. We believe 
that the visible difference between our results and those obtained in Ref.- can be attributed 
to different relativistic effective core potential models used in the calculations. 

Recently long-lived isotopes of the heavier homologue of mercury, Copernicum (Cn or 
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element 112), have been synthesized (see^ and reference therein). At present, the experi- 
ments on Cn are hmited only by the thermochromatography adsorption on the gold surface^; 
several quantum chemical calculations on its compounds and, in particular, Cn2 were per- 
formed to investigate the chemical properties of the new element (including inertness and 
van der Waals radius). The available theoretical data for the Cn dimer (see Table HTl) are 
rather contradictory. The difficulties in the calculations on Cn2 are similar to those for Hg2, 
except for the dramatic increase of the relativistic effects (which, in turn, also affects the 
electron correlations because of the contraction of the s and pi/2 shells). Since there are no 
experimental data on the spectroscopic constants of Cn2, it was important to use the valid 
method tested in calculations on the Hg2 dimer. 

The aim of the present work was to study the role of non-perturbative triple as well as 
of quadruple cluster amplitudes in the bonding of van der Waals dimers Hg2 and Cn2 for 
which accurate calculations of spectroscopic constants are required. (By the effect of non- 
perturbative (or iterative) triples is meant the difference between CCSDT and CCSD(T), 
whereas by that of quadruples is meant the difference between CCSDTQ and CCSDT.) The 
letter is of interest, first of all, from the viewpoint of extracting the van der Waals parameters 
(vdW radius etc.) of the Cn atom. 

II. CALCULATIONS AND DISCUSSION 

The contribution from the iteration of triple and quadruple cluster amplitudes for the 
four valence electrons in Hg2 and Cn2 was evaluated as the difference between the to- 
tal energies obtained in four-electron 4e-FCI and 4e-CCSD(T) calculations^. The ac- 
count of similar contributions from the outer-core electrons is extremely expensive, but 
they are not expected to be noticeable and are not considered here. Scalar relativis- 
tic calculations were performed within the generalized relativistic effective core potential 
(GRECP) model^"— using the MOLCAS^^ code. Generalized correlation bases^*^ com- 
prising the {Us,12p,9d,3f)/[7s,7p,3d,2f] and {16s,2lp,16d,l2f,Ug) /[4:S,6p,3d,2f,lg] 
sets were used for Hg2 and Cn2, respectively. These basis sets were chosen as a com- 
promise between the accuracy and computational efforts in FCI calculations. The ac- 
curacy was checked by calculating the contribution from non-iterative triple cluster am- 
plitudes for the given and rather large {14:S,12p,9d,5f,3g,2h)/[10s,9p,7d,5f,3g,2h] and 



3 



(16s, 21p, 16d, 12 f, lig, lh)/[lls, lOp, 8d, 5/, Ag, Ih] basis sets used in Ref.-^. These calcula- 
tions were performed at the point R = 7.0 a.u. for Hg2 and R = 6.5 a.u. for Cn2 which are 
close to the equilibrium geometry. 

The corresponding contributions from the perturbative (non-iterative) triple cluster am- 
plitudes are —90 cm^^ (large basis set) and —84 cm~^ (small basis set) for Hg2, whereas 
they are —51 cm~^ (large basis set) and —43 cm^^ (small basis set) for Cn2. Then the 
contributions from non-perturbed triples and quadruples for both dimers were evaluated for 
internuclear distances from 5 to 20 a.u. at a step of 1 a.u. and added to the best potential 
energy curves obtained in^ by the CCSD(T) method with large basis sets taking into account 
the counterpoise and spin-orbit corrections. The spectroscopic constants derived from these 
curves are collected in Table [Hand Table [ITl One can see in Table [T] that the agreement with 
the experimental data is significantly improved. Table IIIII presents the energy differences 
E{7.0 a.u.) — E{oo) for Hg2 and £'(6.5 a.u.) — E{oo) for Cn2 (the chosen finite distances are 
close to equilibrium values of the best theoretical level) calculated at the SCF, 4e-CCSD, 
4e-CCSD(T), 4e-CCSDT, and 4e-FCI levels of theory using small basis sets. One can deduce 
from these data that the total contributions from the connected quadruples and from the 
iteration of triples are —23 cm~^ and —6 cm~^ for Hg2 and Cn2, respectively. Since the 
contributions from only the iteration of triples are —20 cm~^ and —5 cm~^, we can conclude 
that the effect of the iterative accounting for the triple cluster amplitudes is noticeable, 
whereas the accounting for the quadruples is not. The potential energy curves for the Hg2 
and Cn2 dimers obtained in SCF calculations are repulsive. They remain repulsive also if 
only valence s shells are correlated. An interesting feature that can be seen from the data in 
Table IIIII is that the potential energy curve of the Cn2 dimer becomes even more repulsive of 
only the valence 7s shell is correlated. To corroborate the reliability of the DFT-based spin- 
orbit corrections from Ref. obtained with the generalized-gradient approximations for the 
exchange-correlation functional, i.e., Perdew-Wang (PW91)^ and Perdew-Burke-Ernzerhof 
(PBE)^^ models, we performed similar calculations with the meta-GGA TPSS functional^. 
All three correcting functions were found to be nearly identical. 
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III. CONCLUSIONS 



The effects of the non-perturbative treatment of triples and the inclusion of quadruples 
on the computed Hg-Hg and Cn-Cn van der Waals interaction energies are evaluated. 
For Hg2, the account of these contributions substantially improves the CCSD(T) results, 
bringing the computed dissociation energies, equilibrium distances, and vibrational constants 
into nearly perfect agreement with the experimental values. One can suppose that the non- 
perturbative treatment of triples can play important role in the accurate description of other 
van der Waals complexes. 
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TABLE I: Calculated equilibrium distance (in A) and other spectroscopic constants (in cm^^) for 
the ground state of the Hg2 molecule as compared to the experimental and previously calculated 
data. 



method 


Re 




We 


WeXe 


DFT(B88/P86)^ 


3.63 


73 


14 




DFT(PW91)^ 


3.55 


385 


24 




'^CCSD(T)3. 


3.60 


581 






CCSD(T) +AEsc^- 


3.729 


379 


19.4 


0.24 


CCSD(T)i 


3.718 


379 


19.4 


0.24 


CCSD(T)i^ 


3.836 


315 


16.3 




CCSD(T) +AEscA 


3.743 


328 


18.4 


0.28 


CCSD(T) +AEs(^ 


3.730 


355 


18.74 


0.24 


This work 


3.711 


377 


19.32 


0.24 


Experiments^ 


3.69 ±0.01 


380 ± 25 


19.6 ±0.3 


0.25 ±0.05 



^ Total relativistic version of the coupled-cluster method. Small basis set, ESSE is not 
compensated. 
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TABLE II: Calculated equilibrium distance (in A) and other spectroscopic constants (in cm ^) 
for the ground state of the Cn2 molecule as compared to the previously calculated data. 



method R^, Df. We We^e 

DFT(B88/P86}^ 3.45 315 25 ??? 

DFT(PW91}^ 3.39 649 30 ??? 

"CCSD(T)3. 3.07 1508 

CCSD(T) +/\Escr- 3.323 768 29.83 0.3272 

This work 3.322 775 29.94 0.3266 



Total relativistic version of the coupled-cluster method. Small basis set, ESSE is not 
compensated. 

TABLE III: Calculated energy differences (in cm^i) E{R = 7.0 a.u.) - E{R = 200.0 a.u.) for Hg2 
and E{R = 6.5 a.u.) - E{R = 200.0 a.u.) for Cn2 



Hg2 Cn2 

SCF 617.4 944.5 

4e-CCSD 211.9 1573.3 

4e-CCSD(T) 128.2 1530.0 

4e-CCSDT 108.3 1524.7 

4e-FCI 105.4 1523.9 
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